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Abstract 



A new evaluation of the constraint on the number of light neutrino species 
(N w ) from big bang nucleosynthesis suggests a discrepancy between the pre- 
dicted light element abundances and those inferred from observations, unless 
the inferred primordial 4 He abundance has been underestimated by 0.014± 
0.004 (la) or less than 10% (95%C.L.) of 3 He survives stellar processing. 
With the quoted systematic errors in the observed abundances and a conser- 
vative chemical evolution parameterization, the best fit to the combined data 
is N„ = 2.1 ± 0.3 (la) and the upper limit is N u < 2.6 (95% C.L.). The data 
are inconsistent with the Standard Model (N v = 3) at the 98.6% C.L. 
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Along with the Hubble expansion and the cosmic microwave background radiation, big 
bang nucleosynthesis (BBN) provides one of the key quantitative tests of the standard big 
bang cosmology. The predicted primordial abundances of 4 He, D, 3 He, and 7 Li have 
been used to constrain the effective number of light neutrino species (N w ) [] The 
neutrino counting includes anything beyond the Standard Model [such as a right-handed 
(sterile) neutrino] that contributes to the energy density. This constraint is complementary 
to neutrino counting from the invisible width of Z decays (Nf), which is sensitive to a much 
larger mass range (^5 Mz/2, where Mz is the Z mass), but only to neutrinos fully coupled 
to the Z; the current result is = 2.988 ±0.023 |J, in agreement with the Standard Model 
(Nf = 3). 

The primordial 4 He abundance is sensitive to the competition between the early universe 
expansion rate and the weak interaction rates responsible for the interconversion of neutrons 
and protons. The expansion rate depends on the overall density and hence on N„, while the 
weak rates are normalized via the neutron lifetime. Recent improvements in neutron lifetime 
measurements have significantly reduced the uncertainty in the 4 He prediction and, coupled 
with increasingly accurate astronomical data on extragalactic 4 He, have led to tighter con- 
straints on N„; at 95% C.L. N„ < 4 in 1989 §, < 3.3 in 1991 §, and < 3.04 in 1994 [§. 
However, a constraint as strong as N„ < 3.04 hints that the standard theory with = 3 
may not provide a good fit to the observations. 

In this Letter we present new BBN limits on N„ and the baryon-to-photon ratio (rj) from 
simultaneous fits to the primordial 4 He, D, 3 He and 7 Li abundances [hereafter we use the 
notation Y p ( 4 He mass fraction), y 2p = D/H, y 3p = 3 He/H, and y 7p = 7 Li/H, fractions by 
number] inferred from the astrophysical observations. In particular, we incorporate new con- 
straints on y 2p 0, which are based on a generic chemical evolution parameterization || and 



Neglecting the baryon contribution, the total energy density ptot depends on N u as ptot = 
Pj + p e + N u p u , where p 7 , p e , and p v are the energy density of photons, electrons and positrons, 
and massless neutrinos (one species), respectively. 
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which significantly improve the prior constraints 0,||J. Our likelihood analysis systematically 
incorporates the theoretical and observational uncertainties. The theoretical uncertainties 
and their correlations are estimated by the Monte Carlo method |n],[|[n],[nj . Non-Gaussian 
uncertainties in the observations, such as the adopted systematic error in the value of Y p , 
the upper and lower limits for D, and the model-dependent 3 He survival parameter (g 3 ), are 
treated in a statistically well-defined way. 

We adopt a primordial helium abundance estimated from low metallicity HII regions 



Y p = 0.232 ± 0.003 (stat) ± 0.005 (syst), (1) 
assuming a Gaussian distribution for the la statistical uncertainty and a flat (top hat) 



distribution with a half width of 0.005 for the systematic uncertainty [12]. The systematic 



error is similar to that used for previous estimates on an d to that obtained from 



Pagel's analysis of the data ||14 |. 



New D constraints were obtained in Refs. |7]||, using pre-solar abundances of D and 3 He 
(as inferred from 3 He measurements in the solar wind, meteorites, and lunar soil and 
a generic chemical evolution parameterization: 

y 2p = (1.5 - 10.0) x 10- 5 (2) 
y 3p < 2.6 x 10~ 5 (95% C.L.). (3) 

Although these constraints are independent of any specific model for primordial nucleosyn- 
thesis, standard BBN or otherwise, they do depend on the adopted 3 He survival fraction 
g%. To be consistent with prior analyses we adopt = 0.25 PJT^Jl[] although the effec- 
tive #3 of most models is significantly larger than this (see later discussion). When the 
observational bounds in Eqs. 2 and 3 are convolved with the BBN predictions (which are a 
function of 7] with N„ fixed at 3), even tighter constraints on D and 3 He may be inferred 
y 2p = (3.5±H) x lO" 5 and y 3p = (1.2 ±0.3) x 10~ 5 at 95% C.L.. The resulting upper bound 
to 7/2 P is roughly 30% lower than the corresponding bound in Ref. Q and this has the effect 



of raising the lower bound on the allowed range of rj. Our central value for y2 P is an order of 
magnitude smaller than the abundance inferred from a possible D detection in absorption 
against a high redshift QSO fIT|]IBfl , but consistent with that reported for a different QSO 



absorption system ||19|| . 

We estimate the primordial 7 Li abundance from the metal-poor stars in our Galaxy's 
halo: 

y 7p = (1.2j£jj) x 10- 10 (95% C.L.). (4) 



This estimate is consistent with other recent determinations pO| , |ilf which take into account 
possible post big bang production and stellar depletion of 7 Li. 

For standard (N u = 3) BBN, the theoretical predictions with the uncertainties (la) 
determined by the Monte Carlo technique are displayed as a function of rj in Fig. [1]. Also 
shown in Fig. |1] are the constraints obtained by our likelihood analysis of the predictions 
and observations. The result is disturbing: the constraints on rj from the observed He and 
D— 3 He abundances appear to be mutually inconsistent. 

To explore this more carefully, all four elements are fit simultaneously, yielding the 
likelihood function for N„ shown in Fig. |] (where the likelihood is maximized with respect 
to rj for each N„). The BBN predictions for the D, 3 He, and 7 Li abundances are sensitive to 
the baryon-to-photon ratio rj, but only weakly dependent on N„. The BBN prediction for 
4 He is very weakly dependent on rj and is approximately proportional to (N^ — 3). In our 
likelihood analysis, we have computed the Monte Carlo predictions for all of the element 
abundances for 1.5 < N u < 4 and 10~ 10 < rj < 10~ 9 . The N u and 77 dependences of the 
uncertainties, the 77 dependence of the correlations among the uncertainties pl|j5|JT^|, and 
the correlations between rj and the y2 P and y% p values have all been included in the likelihood 
function. 

Fig. ^| shows that the Standard Model (N u = 3) yields an extremely poor fit. The best 
fit is for Nj, = 2.1 ± 0.3, and the upper-limit from the joint likelihood (Fig. is 

N, < 2.6 (95% C.L.). (5) 



The ratio of the likelihood of = 3 to the best fit N„ = 2.1 is 0.014. This value provides 
an estimate of the goodness-of-fit of the standard (N u = 3) theory. The result of our 
simultaneous fit in the rj — N u plane is shown in Fig. |^. The constraint on the baryon-photon 
ratio is i] = (4.4ioJ) x 10~ 10 (la). The conflict between the lower and upper bounds on rj 
coming from D and 4 He, respectivley, has been noted before |22|| . Our results exacerbate this 
discrepancy to roughly a 3 standard deviation effect, mainly due to our new D constraint. 

In setting limits when the likelihood function extends beyond the physical parameter 
space, it is usually a reasonable (and conservative) prescription to renormalize the probabil- 
ity density distribution within the physical part of parameter space. This implies that one 
should renormalize the likelihood function for > 3, when constraining any (nonstandard) 
particle contribution in addition to three massless neutrinos in the Standard Model. Exam- 
ining the N„ limit this way, the 95% CL. limit for Nj, extends to 3.25 (for r/ = 4.6 x 10~ 10 ). 
However, we do not advocate this interpretation, since the poorness of the Nj, = 3 fit makes 
this additional constraint for > 3 meaningless. 

The combined data (D, 3 He, 4 He, and 7 Li) with the adopted uncertainties are inconsistent 
with standard (N u = 3) BBN, for a conservative choice of 3 He survival factor g 3 = 0.25. But 
what if some of the uncertainties have been underestimated? In particular, the systematic 
uncertainty in the 4 He observational data may be 3 or more times larger |23| than the 
estimate in Ref. |13[ . With rj determined by the combined D— 3 He and 7 Li constraints, BBN 
predicts Y p = 0.246± 0.002 (la), where the error includes the uncertainties from the D— 3 He 
and 7 Li constraints and from the BBN theory calculation. This value for Y p required for 
BBN consistency is 0.014 above the adopted observed value [Eq. (1)]. 

In Fig. H we show the rj — N u constraints when the central value for Y p is systematically 



2 There is no standard procedure to estimate the goodness-of-fit when non-Gaussian uncertainties 
are involved in a likelihood analysis. In addition to using the ratio of the likelihoods for Nj, = 2.1 
and 3, we have also estimated the goodness-of-fit with the standard x 2 method by approximating 
the errors with Gaussian distributions: the results from the two methods are consistent fll2| . 
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shifted by AY. To be consistent with N„ = 3, AY has to be significantly larger than the 
systematic error adopted in Eq. 1. When AY is fit as a free parameter with N„ fixed to 3, 
we obtain AY = 0.014 ± 0.004 at la. Even allowing AY to change freely, the 7 Li and ISM 
D constraints still bound rj from above at 6.3 x 10~ 10 (95% C.L.); ISM D alone bounds r\ 



from above at 9 x 10 10 . The claim in Ref. [23] that r\ can be as large as ~ 14 x 10 10 is 



unjustified. 

We have also examined (Fig. ]|) how the r]— N u constraint is relaxed when the 3 He survival 
factor, which affects the upper limit on y2 P , differs from that adopted (#3 = 0.25). Relaxing 
the ?/2p upper limit so as to be consistent with the Y constraint requires a significantly 
smaller g 3 . When g 3 is allowed to be a free parameter with N„= 3 fixed, we obtain g 3 < 0.10 
at 95% C.L., i.e. stellar destruction of 3 He would need to be significantly larger than is 
implied by stellar and chemical evolution models. Although it is difficult to assign statistical 
probabilities to various values of g%, one can assess the current status of models of Galactic 
chemical evolution and their associated 3 He destruction. In this Letter we have adopted 
an effective^ g 3 = 0.25, a choice based on the fact that g 3 > 0.25 for any star [|T6"],|9"|,|T| . 
Recent studies ]2j]j2j)| ,[8|] have effective g 3 's larger than 0.25, a fact supported by Ostriker 
and Schramm's analysis of horizontal branch stars which concludes that g% > 0.3 and 
Rood, Bania and Wilson's observation of 3 He in planetary nebulae which suggests that low 
mass stars are net producers of 3 He | 27fl . In order for the effective $3 to be lower than 
0.25, gas would have to be cycled thru several generations of relatively massive stars (which 
are the most efficient destroyers of 3 He) without overproducing metals. Allowing stellar 
3 He production (as evidenced in low mass stars) would effectively increase g 3 and therefore 
exacerbate the present discrepancy between theory and observations. There are models 
and parameterizations which attempt to address these issues. The models of Olive et al. 
include stellar 3 He production in low mass stars and therefore tend towards large g$. 



3 The 53 used in previous BBN analyses is an effective 53 in that it represents the 53 per star 
integrated over all stars and cycled thru some number of stellar generations. 



6 



They conclude that "the only way to reduce g% below that of the massive stars (around 
0.3) would be to argue that the gas in the region has been cycled through stars several 
times. Such an assumption however would invariably predict 4 He abundance factors of 2-4 
higher than those observed." Vangioni-Flam and Casse |29| find that the effective #3 can be 
small, but the associated metals overproduction requires the revision of classical models of 
chemical evolution {e.g., including metal depletion by outflow). The interplay between the 



lower bound to #3 and metal overproduction is reflected in Copi, Schramm and Turner's |3Q 



'stochastic history' parameterization of chemical evolution. Their 95% CL. lower bounds to 
r) are greater than or equal to ours provided they satisfy the metallicity constraint. It is our 
conclusion that our D constraint is robust and probably overly conservative - most models 
of chemical evolution yield D constraints which make the fit between theory and observation 
for Nj, = 3 worse than we report here. For example, if we assume that g^ is equally likely 
to be between 0.25 and 0.5, standard BBN would be ruled out at the 99.1% C.L.. 

The standard (N„ = 3) BBN predictions for the primordial 4 He and D abundances appear 
to be inconsistent with those inferred from observations, unless the inferred primordial 4 He 
mass fraction has been underestimated by AY = 0.014 ± 0.004 or the 3 He survival fraction, 
#3, is smaller than 0.10. While it may be that the crisis lies in the observational data and/or 
its extrapolation to primordial abundances, it is possible to alter standard BBN in order 
to reduce the 4 He prediction to the level consistent with the D constraint. The effective 
N u can be reduced to the range 2.1 ± 0.3 in several ways: massive tau neutrinos, neutrino 
degeneracy, or new decaying particles to name but a few. 

This work is supported by the Department of Energy Contract No. DE-AC02-76-ER01545 
at Ohio State University and DE-AC02-76-ERO-3071 at The University of Pennsylvania. 
R.J.S. is supported in part by NASA (NAG 5-2864). We thank C. Copi, D. Schramm, and 
M. Turner for spirited discussions relating to this paper and others. 
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FIGURES 

FIG. 1. BBN predictions (solid lines) for Y p , y2 P , and yj p with the theoretical uncertainties 
(lcr ) estimated by the Monte Carlo method (dashed lines). Also shown are the regions constrained 
by the observations at 68% and 95% C.L. (shaded regions and dotted lines, respectively). 

FIG. 2. The likelihood function for N„ when the observations for Y p , y2 P , y% p , and y-j p are 
fit simultaneously. For each N„ the likelihood function is maximized for -q. The upper limit is 
N„ < 2.6 (95% C.L.) The fit for the Standard Model (N„ = 3) is excluded at 98.6% C.L. 

FIG. 3. The combined fit of the observations to N u and 7710 = rj 10 10 . 

FIG. 4. The combined fit of the observations when the systematic uncertainty in the 4 He 
observation (AY sys ) is fixed to 0, 0.005, 0.010, and 0.015. 

FIG. 5. The combined fit of the observations when the 3 He survival factor (53) is fixed to 0.10, 
0.25, and 0.50. 
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